A POSTERIORI L°°(i 2 )-ERROR BOUNDS IN FINITE ELEMENT 
APPROXIMATION OF THE WAVE EQUATION 

EMMANUIL H. GEORGOULIS, OMAR LAKKIS, AND CHARALAMBOS MAKRIDAKIS 

Abstract. We address the error control of Galerkin discretization (in space) 
of linear second order hyperbolic problems. More specifically, we derive a pos- 
teriori error bounds in the L°°(L 2 ) -norm for finite element methods for the 
linear wave equation, under minimal regularity assumptions. The theory is 
developed for both the space-discrete case, as well as for an implicit fully dis- 
crete scheme. The derivation of these bounds relies crucially on carefully con- 
structed space- and time-reconstructions of the discrete numerical solutions, 
in conjunction with a technique introduced by Baker (1976, SIAM J. Numer. 
Anal., 13) in the context of a priori error analysis of Galerkin discretization of 
the wave problem in weaker-than-energy spatial norms. 



1. Introduction 

In computing approximate solutions of evolution initial-boundary value prob- 
lems mesh-adaptivity plays an important role, in that it drives variable resolu- 
tion requirements, thereby contributing reduction in computational cost. Adaptive 
strategies are often based on a posteriori error estimates, i.e., computable quantities 
which estimate the error of the finite element method measured in a suitable norm 
(or other functionals of interest). 

A posteriori error bounds are well developed for stationary boundary value prob- 
lems (e.g., |Ver96l lAOOOl iBSOll ICB02l IDM1 ISte07l ICKNS08j and the references 
therein). Adaptivity and error estimation for parabolic problems has also been an 
active area of research for the last two decades (e.g., [EJ951 IVer031 IPic981 HS01, 
IMN031 IBBM051 IBV041 ILM061 IGLV08] and the references therein). 

Surprisingly, there has been considerably less work on the error control of fi- 
nite element methods for second order hyperbolic problems, despite the substantial 
amount of research in the design of finite element methods for the wave problem 
(e.g., |Bak761 IBB791 IBDS79I IBD801 IDS8I1 IJoh93l IMak92l IBJR901 ICJT93I IBJT001 
KM05 and the references therein). A posteriori bounds for standard implicit time- 
stepping finite element approximations to the linear wave equation have been pro- 
posed and analyzed (but only in very specific situations) by Adjerid |Adj02| . Also, 
Bernardi and Siili IBS05] derive rigorous a posteriori L°°(i? 1 )-error bounds, using 
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energy arguments, for the same fully-discrete method we analyze in this work. We 
note that goal-oriented error estimation for wave problems (via duality techniques) 
is also available |BR99[ IBR01] , while some earlier work on a posteriori estimates 
for first order hyperbolic systems have been studied in the time semidiscrete set- 
ting |MN06j . as well as in the fully discrete one |Joh931 [S961IS99] . 

In this work, we derive a posteriori bounds in the L°°(Z/ 2 )-norm of the error, 
which appear to be unavailable in the literature so far. The theory is developed 
for both the space-discrete well as for the practically relevant case of an 

implicit fully discrete scheme. 

The derivation of these bounds relies crucially on reconstruction techniques, used 
earlier for parabolic problems [MN03, LM06, AMN06 . Another key tool in our 
analysis is the special testing procedure due to Baker |Bak76] , who used it in the a 
priori error analysis of Galerkin discretization of the wave problem in weaker-than- 
energy spatial norms. 

While for the proof of a posteriori bounds for the semidiscrete case, the elliptic 
reconstruction previously considered in |MN03l ILM06j suffices, the fully discrete 
analysis necessitates the careful introduction of a novel space-time reconstruction, 
satisfying a crucial local vanishing moment property in time. Our approach is 
based on the one-field formulation of the wave equation and, thus, non-trivial three- 
point time reconstructions are required. A further challenge presented by the wave 
equation is the special treatment of deriving bounds for the "elliptic error" of the 
reconstruction framework, to obtain practically implementable residual estimators. 
The derived a posteriori estimators are formally of optimal order, i.e., of the same 
order as the error on uniform space- and time-meshes. 

The rest of this work is organized as follows. In SJ2]we present the model problem 
and the necessary basic definitions along with the finite element methods for the 
wave equations considered in this work. In <j3] we consider the case of a posteriori 
bounds for the space-discrete problem. In <j4j we derive abstract a posteriori error 
bounds for the fully-discrete implicit finite element method, while in f}5] the case 
of a posteriori bounds of residual type are presented. In ^ we draw some final 
concluding remarks. 



2.1. Model problem and notation. We denote by L p (oj), 1 < p < +oo, to C M. d , 
the Lebesgue spaces, with corresponding norms || • ||lp( w ). The norm of L 2 (uj), 
denoted by ||-|| w , corresponds to the L 2 (w)-inner product (v)w We denote by 
H s (uj), the Hilbertian Sobolev space of order s > of real- valued functions defined 
on to C K d ; in particular Hq (oj) signifies the space of functions in i? 1 (w) that vanish 
on the boundary dco (boundary values are taken in the sense of traces). Negative 
order Sobolev spaces H~ s (u), for s > 0, are defined through duality. In the case 
s = 1, the definition of (•, -) w is extended to the standard duality pairing between 
_ff _1 (o;) and Hq(lj). For 1 < p < +oo, we also define the spaces L p (0,T, X), with 
X being a real Banach space with norm || • consisting of all measurable functions 
v : (0,T) -> X, for which 



2. Preliminaries 



(2.1) 




{ J \\v(t)\\ p x dt) < +oo, for l<p<+oo 
esssup < t < T < +oo, for p = +oo. 



IMIl°°(o,t ; x) : 
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Let f2 C M. d be a bounded open polygonal domain with Lipschitz boundary d£l. 
For brevity, the standard inner product on L 2 (£l) will be denoted by (•, •) and the 
corresponding norm by ||-|| . 

For time t £ (0, T], we consider the linear second order hyperbolic initial- 
boundary value problem of finding u £ L 2 (0, T; Hq(Q)), with u t £ L 2 (0, T; L 2 (Vt)) 
and u tt G I?{f),T;H- x {tl)) such that 

(2.2) u tt - V • (aVti) = / in (0, T) x ft, 

where / G L 2 {Q,T;L 2 {Q)) and a is a scalar-value function in G C(f2), with < 

ttmin < a < «max, Such that 

u(x, 0) = u (x) on ft x {0}, 

(2.3) u t (x, 0) = ui(x) on fix {0} 

u(0,t) = 0on9fix (0,T], 

where u G #o(ft) and Ui G L 2 (fl). 

We identify a function v £ ft x [0, T] -> K with the function u : [0, T] -> (fi) 
and we use the shorthand v(t) to indicate v(-,t). 

2.2. Finite element method. Let T be a shape-regular subdivision of ft into dis- 
joint open simplicial or quadrilateral elements. Each element k £ T is constructed 
via mappings F K : k — ?> k, where k is the reference simplex or reference square, so 
that = U ksT k [Cia78j . 

For a nonnegative integer p, we denote by V v {k) either the set of all polynomials 
on k of degree p or less, when k is the simplex, or the set of polynomials of at most 
degree p in each variable, when k is the reference square (or cube). We consider p 
fixed and use the finite element space 

(2.4) V h := {v £ flJtfi) : v\ K o F K £ V p (k), K £ T}. 

Further, we denote by T := U K , e j-(dK\dft), i.e., the union of all (d— l)-dimensional 
element edges (or faces) e in fi associated with the subdivision T excluding the 
boundary. We introduce the mesh-size function h : f2 — > R, defined by /i(x) = 
diamK, if x £ k and h(x) — diam(e), if x £ e when e is an edge. 

The semidiscrete finite element method for the initial-boundary value problem 
(|23] ) -(|23 ]) consists in finding U £ L 2 (0, T\ V h ) such that 

(2.5) (U„,V) + a(U,V) = (f,V) W G L 2 (0, T; Vh), 
where the bilinear form a is defined for each z,v £ Hq(O) by 

(2.6) a(z,v) — / aVz-Vvdx, 

and the corresponding energy norm is defined for v £ Hq(CI) by 

(2.7) |H| = ||V5V«||. 

To introduce the fully-discrete implicit scheme approximating (|2.2I) - (|2.3[) . we 
consider a subdivision of the time interval (0,T] into subintervals (t n ~ 1 ,t n ], n = 
1, . . . , TV, with t° = and t N = T, and we define k n := t n — t n ~ x , the local time-step. 
Associated with the time-subdivision, let T", n = 0, . . . , N, be a sequence of meshes 
which are assumed to be compatible, in the sense that for any two consecutive 
meshes 7 - ™" 1 and T™, T™ can be obtained from 7 - ™ -1 by locally coarsening some 
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of its elements and then locally refining some (possibly other) elements. The finite 
element space corresponding to T™ will be denoted by V£ ■ 

We consider the fully discrete scheme for the wave problem (|2.2[) . (|2.3p 

for each n = 1, . . . , N, find U n £ V£ such that 

(d 2 U n , V) + a{U n , V) = (/", V) W £ V£, 

where f n := f(t n ,-), the backward second and first finite differences 

dU n - dU 11 - 1 



(2.8) 



(2.9) d z U n := 

with 



jjn _ jjn-1 



for n= 1,2, ...,N, 



(2.10) dU n := { k n 

_V :=ir° Ul for ?i = 0, 

where U :— ir°uo, and tt° : L 2 (il) —> V® a suitable projection onto the finite 
element space (e.g., the orthogonal L 2 -projection operator). 

3. A POSTERIORI ERROR BOUNDS FOR THE SEMI-DISCRETE PROBLEM 

We derive here a posteriori error bound for the error \\u — J7||L°°(o.T;L 2 (a)) be- 
tween the exact solution of (|2.2p . (|2.3[) and that of the scmidiscrete scheme l2~5l 



3.1. Definition (elliptic reconstruction and error splitting). Let U be the (semidis- 
crete) finite element solution to the problem (|2.5|) . Let also II : L 2 (fl) — > Vh be the 
orthogonal L 2 -projection operator onto the finite element space Vh- We define the 
elliptic reconstruction w — w(t) £ Hq(Q), t £ [0, T], of U to be the solution of the 
elliptic problem 

(3.1) a(w,v)={g,v) Vt> G f#(fi) 
where 

(3.2) g := AU-U/ + /, 

and A : Vh —> Vh is the discrete elliptic operator defined by 

(3.3) for q€V h , (Aq, X ) = a(q, X ) V X 6 V h . 
We decompose the error as follows: 

(3.4) e U — u = p — e, where e := w — U, and p :— w — u. 

3.2. Lemma (error relation). With reference to the notation in \3.J$ we have 

(3.5) (e tt ,v)+a(p,v)=0 € H^(Q). 

Proof. We have, respectively, 

(e tt ,v) + a{p, v) = (U tt ,v) + a(w, v) - (/, v) 

(3.6) ={U a ,Uv) + a(w,v)-(f,v) 

= -a(U, Tlv) + a(w, v) + (II f -/,«) = 0, 

observing the identity a(U, Hv) — (ILf — /, v) = a(w, v), due to the construction of 
w. □ 
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3.3. Theorem (abstract semidiscrete error bound). With the notation introduced 
in j3.4\), the following error bound holds: 



l|e||L~(o,T;^(n)) <\\e\\ L ~ (0iT . LH n)) + V2(\\u - U(0)\\ + ||e(0)||) 

(3.7) 

Iktll +c a , T \\u l -Ut(o)\\ > 



where G a ,T '■= min{2T, y/2Cn/oi in in}, where Cq is the constant of the Poincare- 
Friedrichs inequality \\v\\ 2 < Cn||Vu|| 2 , for v € Hq(H). 

Proof. We use a testing procedure due to Baker |Bak76j . Let v : [0,T]xfl->l 
with 



(3.8) v(t,- 



J p(s,-)ds, te [0,71, 



from some fixed r € [0, T]. Clearly v £ i7g(f2) as p E Hq(£1). Also, we observe 
that: 

(3.9) v(t,-)=0, Vu(t,-)=0, and v t (t, •) = -p(t, •), a.e. in[0,T]. 

Set v = v in (|3.5[) . integrate between and r with respect to the variable t and 
integrate by parts the first term on the left-hand side, to obtain 



(3.10) - / (e t ,v t ) + (e t (r), fi(r)) - (e t (0), «(0)) + / a(p, 5) - 0. 

Jo Jo 

Using (|3.9p . we have 

(3 - n) r ^ ||p|12 ~ £ \jt a ^ ] = j T{eup) + < et (°)<^' 

which implies 

(3-12) \\\p{r)\\ 2 \\\ P m\ 2 + \a{v{0),m) = ^V?> + M0),fi(0)>. 

Hence, we deduce 
(3.13) 

^IIp(t)|| 2 - \\WW + -0(5(0), S(0)) < max^llpWH jf ||ej + ||e,(0)||||S(0)||. 

Now, we select r = f such that ||p(t)|| = maxo<t<T ||p(^)||: and we present two 
alternative, but complementary, ways to complete the proof. 

In the first way, we start by observing that ||u(0)|| < r||/5(f)|| , gives 

(3-14) \\\p{r)\\ 2 - \\\ P m 2 < (f \\dt4 +r|M0)||) 2 . 

Using the bound ||p(0)|| < ||e(0)|| + ||e(0)|| , e(0) = U(0)-u ande t (0) = Ut(0)-u u 
we conclude that 

||e|U- (0lT .£ 3( n)) <||e|Uoo ( o lT .i» ( n)) + V2(\\u - U(0)\\ + ||e(0)||) 
(3.15) t 

+ l(J o INI +T\\ Ul -dU{0)\\y 

The second alternative, described next, consists in a different treatment of the 
last term on the right-hand side of (|3.13|) . The Poincare-Friedrichs inequality and 
the positivity of the diffusion coefficient a imply |jw(0)|| 2 < Cf2a nl f n ||C(0)|| 2 , for 
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some constant Cn depending on the domain Q only. Combining this bound with 
(|3.13|) . we arrive to 



(3.16) - \\p(r)\\ 2 -- ||„(0)|| 2 <max T ||^)|| ^ ||e t || + -C a c£j e t (0)\\ 
which implies 

+ V2(||«o- 17(0)|| +||e(0)|| 

(3-17) 

|et|| +A/2Cn/a min ||ui-tf t (0)||. 



Taking the minimum of the bounds (|3 . 1 5[) and (13. 17)) yields the result. □ 

3.4. Remark (short and long integration times). The use of two alternative argu- 
ments in the last step of the proof of Lemma l3~2l improves the "reliability constant" 
C a ,T that works for both the short-time and the long-time integration regimes. 

3.5. Remark (Completing the a posteriori estimation). To obtain a practical a 
posteriori bound, we need to estimate the norms involving the elliptic error e. 
By construction, the elliptic reconstruction w is the exact solution to the elliptic 
boundary-value problem (|3.1|) whose hnite element solution is U . Indeed, inserting 
v = V G Vh in (13. ip . we have 

(3.18) a(w, V) = (AU - ILf + f, V) = a(U, V), 

which implies the Galerkin orthogonality property a(w — U, V) = 0. Therefore, 
by construction, e is the error of the finite element method on Vh for the elliptic 
problem 

(3.19) - V ■ (aVw) = ,g, 

with homogeneous Dirichlet boundary conditions, with g defined by 



3.6. Definition. For every element face e C T, we define the jump across e of a 
field w, defined in an open neighborhood of e, by 



(3.20) M(x) = lim (w(x + 5n e ) - w(x - 6n e ) • n e , 

for x € e, where n e denotes one of the two normal vectors to e (the definition of 
jump is independent of the choice). 

3.7. Theorem (elliptic a posteriori residual bounds [Vcr96, AO00 ). Let z £ Hq(H) 
be the solution to the elliptic problem: 

(3.21) - V ■ (aVz) = r 

r G L 2 (Q) and fl convex, and let Z € Vh be the finite element approximation of z 
satisfying 

(3.22) a(Z, V) = (r, V) W € V h . 

Then, there exists a positive constant C c \, independent ofT, h, z and Z , so that 

(3.23) \\z-Z\\ 2 <C cl £(Z,r,T), 
where 

(3.24) ^^^^(^(ll^r + V-CaV^II^+^II^IaVZ]^)) 172 . 

KST ecr 
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3.8. Corollary (scmidiscrcte residual-type a posteriori error bound). Assume that 
the hypotheses of Theorems \3.3\ and \3.7\ hold. Assume further that f is differ entiable 
with respect to time. Then the following error bound holds: 

ll e llL~(o.T;L2(n)) <C c i||£ {U,g,T)\\ L ~( ,T) + 2C c i / £(U t ,g t ,T) 

Jo 

(3 - 25) +V2C cl £(U(0),g(0),T) 

+ V2||uo- 17(Q)|| +C a , T }\u 1 -Ut(0)\\. 
Proof. Using (|3.18l) . ||e|| and ||et|| can be bounded from above using (|3.23l) . □ 

3.9. Remark. A bound of the form (|3.23[) is only required to to hold for Corollary 
13.81 to be valid. Therefore, other available a posteriori bounds for elliptic prob- 
lems |Ver961 lAOOOj can be also used. 

4. A POSTERIORI ERROR BOUNDS FOR THE FULLY DISCRETE PROBLEM 

The analysis of E|3] is now extended to the case of a fully-discrete implicit scheme 
with the aid of a novel three point space-time reconstruction, satisfying a crucial 
vanishing moment property in the time variable. 

4.1. Definition (space-time reconstruction). Let U n , n = 0, . . . , N, be the fully dis- 
crete solution computed by the method (|2.8j) . II" : L 2 {VL) — > Vfi be the orthogonal 
I/ 2 -projection, and A n : — > Vu to be the discrete operator defined by 

(4.1) forqeV?, (A n q, x ) = a(q, X ) V X G V£. 

We define the elliptic reconstruction w n G i2g(Q), of U n to be the solution of the 
elliptic problem 

(4.2) a(w n ,v) = {g n ,v) Vt> G Hq(Q,), 
with 

(4.3) g n := A n U n - II"/" + /", 

where /"(•) := /(0, ■) and /"(■) := k' 1 f(t, -)dt fern = 1, . . . , N. Finally, we 
need to define the elliptic reconstruction dw° G Hq(Q,), of V° to be the solution of 
the elliptic problem 

(4.4) a(dw°,v) = (dg°,v) Vv G £#(fi), 
with 

(4.5) 9 9 °:=iV-nY + / . 

The time-reconstruction U : [0, T] x £1 — > R of {U n }^ =Ql is defined by 

(4.6) U(t) := t —^—U n + 1 -r^U n - 1 - {t - ~ t)2 d 2 U n , 

for t G (t n_1 ,t n ], n= l,...,N, noting that dll° is well defined. We note that U 
is a C 1 -function in the time variable, with U(t n ) = U n and Ut(t n ) = dU n for , 
n = 0,l,...,N. 

We shall also use the time-continuous elliptic reconstruction w, defined by 

(4.7) w(t) := t —l—w n + 1 ~ r ^w n - 1 - ( t~ tn ~ 1 K tn ~ t)2 &w n , 
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noting that dw° is well defined. By construction, this is also a C 1 -function in the 
time variable. 

We decompose the error as follows: 

(4.8) e := U — u = p — e, where e := w — U, and p :— w — u. 

4.2. Remark (notation overload). In this section we use symbols, e.g., U, w, e, e, p, 
that where used in but with a slightly different meaning. Indeed, these are 
now fully-discrete constructs, corresponding in aim and meaning, but different, to 
their semidiscrete counterpart. It is hoped that this overload of notation should 
not create any confusion. 



4.3. Proposition (fully-discrete error relation). For t e (t n ,t n ], n = 1,...,N 
we have 



(4.9) (e tt ,v)+a{p,v) = ((I-n n )U tt ,v)+^(d 2 U n ,Il n v)+a(w~w n ,v) + (f n -f,v}, 

for all v € Hq(Q), with II™ : L 2 (il) — > Vu denoting the orthogonal L 2 -projection 
operator onto V£ , I is the identity mapping in L 2 (Q), and 

(4.10) ^(t)~-6k-\t-t n -i), 
where i n ~5 : = 

Proof. Noting that U tt {t) = (1 + fi n (t))d 2 U n , for t e {t n ~ l ,t n ], n = 1, . . . , N, and 
the identity a{U n ,U n v) - (U n f n - f n ,v) = a(w n ,v), we deduce 

(e u ,v) +a{p,v) = (U tt ,v) +a{w,v) - (f,v), 

= ((I - U n )U tt ,v) + (Uu, U n v) + a(w, v) - (/, v), 

(4.11) = ((i -n n )u tt7 v) + f i n (t)(d 2 u n ,u n v) 

- a(U n , Tl n v) + a{w, v) + (IP f n - /, v) 
= ((I-U n )U tu v)+ fi n (t)(d 2 U n ,n n v) + a (w-w n ,v) + (/" -f,v), 

□ 

4.4. Remark (vanishing moment property). The particular form of the remainder 
p n (t) satisfies the vanishing moment property 

(4.12) / fx n (t) dt = 0, 

which appears to be of crucial importance for the optimality of the a posteriori 
bounds presented below. 

4.5. Definition (a posteriori error indicators). We define in a list form the error 
indicators which will form error estimator the fully discrete bounds. 

mesh change indicator: r/i(r) := 771,1(7") + i]i,2{i~), with 

m-l j-T 

(4.13) mA (r):=J2 -V>)U t \\ + \\ (I U m )U t \\ , 
and 

m— 1 

(4.14) m . 2 (r) := £ (r - **)|| (IP +1 - W)dW\\ + r\\ (I - n°)V°(0)|| , 
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evolution error indicator: 



(4.15) m 



(r) ■= f \\G\\, 
Jo 



where Q : (0,T] -> R with G\( t j-\tJ] ■= Q 3 3 = h • ■ ■ , N and 

fc 2 fc 3 . 

with g J as in Definition 14. 1 1 and 7j := 'fj-i + -^-dg 3 + -^d g 3 , j = 1, . . . ,N, 
with 70 = 0; 
data error indicator: 



(4-17) m (r) := JL £ ( jT W - /l| 2 ) 1/2 + ( jf ^ *4llf" - /II 2 



1/2 / /• r _ > 1/2 



time reconstruction error indicator: 



(4.18) n 4 (r) := ^ £ ( jT W^W) ^ + ( £ ^ ^IIm"^'! 2 " ^ 

4.6. Theorem (abstract fully-discrete error bound). Recalling the notation of Def- 
inition \4-l\ and the indicators of Definition \j.5\ we have the bound 



||e|Uoo(o,tW;£2(n)) <l|e||L~(o,t«;L 2 (^)) + \/2M|uo - 17(Q)|| + ||e(0)|M 
(4.19) r t» 4 

+ 2 (/ N| + E 7 ^)) + c a,jvii«i - V°\\ , 



where C a ,N '■— niin{2t , -\/2Cn/a m m}, Cn is Poincare-Friedrichs inequality con- 
stant. 



Proof. The proof of Theorem l4.6[ is spread in this and the following paragraphs up 
to 

Next we set v — v in (|4.9j) with v defined by (|3.8|) where p is defined as in (|4.8|) 
(i.e., the fully discrete p), assuming that t" 1 ^ 1 < t < t m for some integer m with 
1 < ti < iV. We integrate the resulting equation with respect to t between and 
r, to arrive to 



(4.20) / (e tt ,v)+ a{p,v)=I l {T)+I 2 {T)+X 3 {r)+X A {T), 

Jo Jo 
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where 



(4.21) 



m-l . t i pr 

ii(r)~J2 {{i-w)u tu v)+ ((i-n m )u tt ,v), 

m-l „t j pr 

1*i{f) '■= J a{w — w^v) -\- J a(w — w m , v) 



3=1 

m— 1 .(3 



J - <.t J pT 

Mr) ■= £ / (/' -/.«>+/ <f ~ /. «) . 

m-l „fJ „ T 

Mr)-=y fi j (d 2 u j ,Wv)+ fi m (d 2 u m ,ir 



In Lemmas 14. 71 14.81 14.91 and 14.111 we will derive bounds of the form 

(4.22) Jk(T)<TH(T)mm\\p(t)\U 

for i = 1, 2, 3, 4. With the help of these, we will conclude the proof in §4.121 □ 

4.7. Lemma (mesh change error estimate). Under the assumptions of Theorem \4-6\ 
and with the notation ft4-21ty we have 

(4.23) I 1 (T)<m(r) (j n^ r l|p(*)||. 

Proof. Observing that the projections IP, j — 1,...,N, commute with time- 
differentiation, we integrate by parts with respect to t, arriving to 

m-l . t i . T 
Mr) = J2 ({I-TP)U t ,p)+ ((I-n m )U tlP ) 

(4.24) 3 

y ' m-l 

+ £((if +1 -w)u t (nv(t j )) - ((/-n°)Ui(o),«(o)). 

The first two terms on the right-hand side of (14.241) are bounded by 

/m-l „t j p t 

(4.25) max JlpWII £ / ||(I - IP)[/ t || + / ||(/-ir™)[/ t 

Recalling the definition of v and that U (P) = dU^ , j = 0, 1, . . . , N, we can bound 
the last two terms on the right-hand side of (|4.24p by 



(4.26) max ||p(f)|| £(r - i')||(IF + 1 - IF )5EP"|| + r||(/ - II»)|| 

V 3=1 

□ 

4.8. Lemma (evolution error bound). Under the assumptions of Theorem \4-6\ and 

with the notation we have 

(4.27) l 2 (T)<m(T)max T \\p(t)\\. 
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Proof. First, we observe the identity 

(4.28) w - w j = -{P - t)dw J + (kj l {P - t) 3 - (P - tf\d 2 w> , 
on each (P" 1 , P], j — 2, . . . , m. Hence, from Definition 14. 1 L we deduce 

(4.29) a(w - w j ,v) = {-{P - t)dg J + (kj\P - t) 3 - {P - t) 2 )d 2 g j ,v) 

The integral of the first component in the inner product on the right-hand side 
of (|4.29p with respect to t between (t J ' -1 ,t 5 '] is then given by Q. The choice of 
constants in Q implies that Q is continuous on P , j — 1, 2, . . . , N and Q(0) = 0. 

Hence, integrating by parts on each interval (P , P], j = 1, . . . , m, we obtain 

(4-30) 1 2 (t) = f (Q,p), 

Jo 

which already implies the result. □ 

4.9. Lemma (data approximation error bound). Under the assumptions of Theo- 
rem \4-6\ and with the notation tJJMy we have 

(4.31) l 3 (T)< V3 (r)max T \\p(t)\\. 
Proof. We begin by observing that 

(4.32) f CP - /)=0, 

Jv- 1 

for all j = 1, . . . , m — 1. Hence, we have 

m-l „t j m-l pt j 

(4.33) / (f "/»«>= E / & ~ f> % ~ 

where v 3 { ) := kj J t j-i v(t, -)dt. Using the inequality 

P j k 2 r tj 

(4.34) J \\Z-&f <_L_ jf |M| 2 , 

and recalling that v t = p, we have, respectively, 

m— 1 m—i r t j . i in , /■{ 



1 ^ / \ 1/2 / ft j 

(4.35) <ig (X-. "'■') (/,-,*> 

i '^z 1 / _ \ 1/2 

^E (LW-Sf) -xJIpWII 



27T ^— ' V / H -l 3 J 0<t<T 

For the remaining term in I3, we first observe that 

(4-36) / \\v\\ 2 dt< f k m ! ||p|| 2 d S d<<fc^ max ||p(<)|| 2 , 

which implies 

1/2 



(4-37) / {f m -f,v)<( kUr-fW 2 ) max ||p(t)| 

l+m-l \ lfm-1 I 0<t<T 
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Recalling 773 from Definition 14.51 we conclude the proof. □ 

4.10. Remark (the order of the data approximation indicator). The choice of the 
particular combination of functions involving the right-hand side data / in the 
definition of g n in the elliptic reconstruction, results to the property (|4.32j) . When 
/ is differentiable, we have 773 (r) = 0{k 2 ) as k :— maxi<j< m kj — > 0, and the 
convergence is of second order with respect to the maximum time-step. In this 
case, 773 is, therefore, a higher order term. 

4.11. Lemma (time-reconstruction error bound). Under the assumptions of Theo- 
rem \4-6\ and with the notation l{4.21\ ) we have 

(4.38) Zk(r)<Tfc(T) nu« \\p(t)\\. 

Proof. The method of bounding Z4 (t) is similar to that of Lemma 14.91 so we shall 
only highlight the differences. 

Recalling the vanishing moment property (|4.12p and noting that d 2 U^ is piece- 
wise constant in time, we have 

m-l „t j to— 1 „t j 

(4.39) n j {d 2 U j ,n j v)=y2 if {&TP,W 
j= i Jt3 - 1 3=1 JtJ - x 

— ■ 1 r*t^ 

where v J (■) — k, Jp-i v(t, -)dt. Hence, since IP commutes with time integration, 



we obtain 



(4.40) 

TO-l -t^ , m-l „ t j 1/2 . rt j 

Y, / u?{d 2 w,w{v-P)) < — J2 ( / IIm^ 2 ^'II 2 ) ( /. fc 2 l|n-V|| 2 ) 

. m — 1 „t j 

^Eiy^^ii^ii 2 ) ^ r iip(*)n. 

For the remaining term in Z 4 , upon using an argument similar to (I4.36p . we have 

(4.41) J T m Jn m d 2 U m ,U m v) < (jT_ t eilM m 5 2 f/ m |r) 1/2 o max ||p(i)||. 

Recalling the definition of 774 in ^4.51 we conclude. □ 

4.12. Concluding the proof of Theorem 14.61 Starting from (|4.20[) . integrating 
by parts the first term on the left-hand side, and using the properties of v, we arrive 
to 



1/2 



which implies 



r 4 

J o i=l 



(4-43) I||p( T )|| 2 -i|| /5 (0)|| 2 + i a (5(0),«(0)) = 



r 4 

Jo 
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Hence, we deduce 

(4.44) i||p(r)|| 2 -i||p(Q)|| a + ia(t;(0),t5(0)) 



< max \\p(t)\\([ |H| + £>(r) 
— Jo i=i 



|e t (0)||||5(0)|| 



We select r = f such that ||/c(f)|| = max < t < t iv ||p(i)|| . First, observing that 
||5(0)|| <r||p(f)||, gives 

(4-45) \\\p{rW \\\pm\ 2 < ( f Ml + £>(r) + r||e t (0)|| 

~' t=i 

Using the bound ||p(0)|| < ||e(0)|| + ||e(0)|| and observing that e(0) = U(0) - tt(0) 
U° — uq and that e*(0) = Ut(0) — ut{0) = V a — u\, we arrive to 



e||L-(o,t«;L 2 (o)) <l|e||L~(o,t«;L 2 (o)) + V^H^o - *7°|| + ||e(0)||J 

+ 2 ( / IMI+J>(t 



(4.46) . ,t N 4 

' ' " < A ' 1 t N \\ Ul -V° 

The second way is completely analogous to the proof of the semidiscrete case. 
5. Fully-discrete a posteriori estimates of residual type 

To arrive to a practical a posteriori bound for the fully-discrete scheme from 
Theorem 14.61 the quantities involving the elliptic error e should be estimated in 
an a posteriori fashion: this is the content of Lemmas 15.11 and 15.31 below, when 
residual- type a posteriori estimates are used. 

5.1. Lemma (estimation of the elliptic error). With the notation introduced in 
Definition ^. 1\ we have 

(5.1) INU~(o,t";Z'(n)) + V2||e(0)|| < S x (f) + V2C el £°, 

where 



S 1 (t N ) :=max{^C cl £0/°,%°,T o ), 



(5.2) 



35 31 kj 

— max - 



(C*£>' + Caa^Jf* - f\\)}, 



,27 27 i<j<N kj-i ) o<j<N 

with & := £{W,AW j - IP p + f j ,T j ), j = 0,l,...,N. 

Proof. For t £ (t j ~ l ,t' J ], j = 1, . . . , N, we have 
(5-3) 

e = —T——( wj - U j ) + t -^(w 1 - 1 - W- 1 ) - {t ~ P 1)(P ' f)2 (9V - d 2 W), 
kj kj kj 

from which, we can deduce 

(5.4) ||e|| <maxf(^ + ^ max -^-) max \\w> - W || , ^\\dw a - V°\\ }, 
y ' 11 11 - l\27 27 i<j<N kj-\J o<j<N n '"27" "J 

noting that 

{t-t^){V-tf 4fcf 

(5.5) max = — -. 

t£{ti- l ,ti] kj 27 
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It remains to estimate the terms |] w J — U J \\ and \\dw° — V°\\. To this end, recalling 
the notation of Definition l4.il we define wl E Hq(Q) to be the solution of the elliptic 
problem 

(5.6) a(wi,v) = (AU 3 - IP f 3 + f 3 \v) Vv e ftf 

for j = 0, 1, . . . , N. Note that, due to the fact that /° = /°, we have = w°. 
By construction, we have a(wt,V) = (A 3 U 3 - U> f 3 + f 3 ,V) = a(U 3 ,V) for all 
V G Vjl , j = 0, 1, . . . , N. Hence, U 3 is the finite element solution (in V£) of the 
elliptic boundary- value problem (|5.6[) . In view of Theorem 13 .71 this implies that 



(5.7) \\wl-U 3 \\<C cl £ 3 , 

for j = 0, . . . , N . Similarly, by construction, we have a(dw°, V) = (A°V° — II / + 
/°, V) = a{V°, V) for all V € Vg. Hence, 

(5.8) \\dw° - dU°\\ < C e i£(V°, dg°, T°). 

Moreover, since w 3 — wl is the solution of an elliptic problem with right hand-side 
f 3 — f 3 , standard elliptic stability results yield 

(5-9) \\w 3 -wi\\ <Cna-U\f 3 -f 3 l 

for j = 1, . . . , N. Finally, using the triangle inequality 

(5.10) \\w 3 - U 3 \\ < \\w 3 -wi\\ + \\wt - U 3 \\, 

along with the bounds (|5.9[) . (|5.8|) and (|5.7I) . already implies the result. □ 



5.2. Remark. The bound (|5.1[) contains both the elliptic estimators £ (•, •, •) and the 
data-oscillation terms \\f 3 — f 3 \\ which are, in general, of first order with respect to 
the time-step. The data-oscillation terms are expected to dominate the data error 
indicator 773 (cf. Remark l4.10p . On the other hand, if the numerical scheme (|2.8|) is 
altered so that f 3 — f 3 (as done, e.g., in [Bak76 ), then the data-oscillation terms 
in (|5.ip vanish. Similar remarks apply to the result of Lemma 4.12 below. 

For each n = 1, . . . , N, we denote by T™ the finest common coarsening of T n 
and T™ -1 , and by Vfi :— V£ n V£~ , the corresponding finite element space, along 
with the orthogonal L 2 -projection operator LT" : L 2 (il) — > V£ . 

5.3. Lemma (estimation of the time derivative of the elliptic error). With the 
notation introduced in < j^. 1\ we have 

t N 

(5.11) / |H| < s 2 (t N ), 

where 

2 N / 

(5.12) 6 2 (t N ) := - ^(2% + k i+1 )(C A e> B + C n a m lJdf 3 - df\ 

j=o 

with 

(5.13) £i~S(dW,d(AW)-d(Wf) + df,P), j =0,1,..., TV. 
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Proof. For t G (P j = 1, . . . , N, we have 

(5.14) e t = dw 3 - dW ~ kj l (V - t)(t j - 2t j ~ 1 + t)(d 2 w 3 - d 2 W), 
from which, we deduce 

ft 3 Ah. Oh 

(5.15) / ||e t || < -^\\duP -dU j \\ +^\\dwi- 1 -811*-% 
Jti-i 3 3 

noting that 

f^k-\V-t){V-2t^ + t)= 2 -±L. 



(5.16) 



Combining (15. 15ft for j = 1, . . . , iV, we arrive to 

t N N 

(5.17) / \\e t \\<lj2(2k j+ k j+1 )\\dwl-dWl 

with fco = and kjy+i = 0. 

It remains to estimate the terms \\dw J — 9EP||. To this end, recalling the def- 
inition of the functions wi G Hq(SY) from the proof of Lemma 15.11 and, since 
Vl := V( n V^~\ we have a(wi, V) = a(W , V) for all V G V 3 and a(w^\V) = 
aiUi- 1 , V) for all V G V^' , for j = 1, . . . ,N. Therefore, we deduce 

(5.18) a{dwi, V) = a{dU°, V) for all V G V£, 

for j = 1, . . . , JV, i.e., dU 3 is the finite element solution in Vl of the boundary- value 
problem 

(5.19) a(dwi,V) = (d{A j U j ) - d(W f) + df , v) Wv G J#(fi). 
In view of Theorem 13.71 this implies that 

(5.20) \\dw{-dW\\ <C e i4, 

for j = 1, ,.,,N. We also recall that, by construction, we have a(dw°,V) = 
a(V°, V) for all V G V® . Hence, dEHJ) also holds. 
Moreover, since 

(5.21) a{dw 3 \ V) = (d(A j U j ) - 9(IF f) + df j ,v) Vv e H%(Sl), 
j = 1, . . . , N, (cf. Definition 14. 1| . As in (|5.9|) . elliptic stability implies 

(5.22) \\dw 3 - dwi\\ < Cua^Jdf - df\\ , 
for j = 1, . . . , N and, using the triangle inequality 

(5.23) \\dw i - dW\\ < \\dw 3 - dwi\\ + \\dw{ - dU j \\, 

along with the bounds ()5.22|) . (|5.8|) and (|5.20p . already implies the result. □ 

5.4. Theorem (fully-discrete residual- type a posteriori bound). With the same 
hypotheses and notation as in Theorems \4-6\ and \3.7\ we have the bound 

l|e||L~(o, t « ; ^(o)) <Si(t N ) + V2C cl £° + V2\\u - U(0)\\ 

(5 24) 4 

+ 2S 2 (t N ) +2j2vi(t N ) + C a , N \\ Ul -V°\\, 

i=l 

where Si,£° are defined in Lemma l5.il S2 is defined in Lemma \ 5.3l and r\i, i = 
1,2,3,4 after (41) respectively. 
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Proof. Combining Theorem 14.61 with the bounds derived for e in Lemma 15.11 and 
et in Lemma 15. 3[ we arrive to an a posteriori error bound. □ 

6. Final remarks 

The design and implementation of adaptive algorithms for the wave equation 
based or rigorous a posteriori error estimators is a largely unexplored subject, de- 
spite the importance of these problems in the modeling of a number of physical phe- 
nomena. To this end, this work presents rigorous a posteriori error bounds in the 
L°°(Z/ 2 )-norm £ Qr seconc j order linear hyperbolic initial/boundary value problems. 
The derived bounds are formally of optimal order. The numerical implementation 
of the proposed bounds in the context of adaptive algorithm design for second order 
hyperbolic problems remains a challenge that deserves special attention and will be 
considered elsewhere. 



References 



[Adj02] Slimane Adjerid, A posteriori finite element error estimation for second-order hyper- 
bolic problems, Comput. Methods Appl. Mech. Engrg. 191 (2002), no. 41-42, 4699- 
4719. MR MR1929627 (2003g:65113) 

[AMN06] Georgios Akrivis, Charalambos Makridakis, and Ricardo H. Nochetto, A posteriori 
error estimates for the Crank- Nicolson method for parabolic equations, Math. Comp. 
75 (2006), no. 254, 511-531 (electronic). MR MR2196979 (2007a:65114) 

[AO00] Mark Ainsworth and J. Tinsley Oden, A posteriori error estimation in finite element 
analysis, Pure and Applied Mathematics (New York), Wiley-Interscience [John Wiley 
& Sons], New York, 2000. MR MR1885308 (2003b:65001) 

[Bak76] Garth A. Baker, Error estimates for finite element methods for second order hyperbolic 
equations, SIAM J. Numer. Anal. 13 (1976), no. 4, 564-576. MR MR0423836 (54 
#11810) 

[BB79] Garth A. Baker and James H. Bramble, Semidiscrete and single step fully discrete 

approximations for second order hyperbolic equations, RAIRO Anal. Numer. 13 (1979), 

no. 2, 75-100. MR MR533876 (80f:65115) 
[BBM05] A. Bergam, C. Bernardi, and Z. Mghazli, A posteriori analysis of the finite element 

discretization of some parabolic equations, Math. Comp. 74 (2005), no. 251, 1117—1138 

(electronic). MR MR2136996 (2007c:65072) 
[BD80] Garth A. Baker and Vassilios A. Dougalis, On the L°° -convergence of Galerkin ap- 
proximations for second-order hyperbolic equations, Math. Comp. 34 (1980), no. 150, 

401-424. MR MR559193 (81f:65066) 
[BDS79] Garth A. Baker, Vassilios A. Dougalis, and Steven M. Serbin, High order accurate two- 
step approximations for hyperbolic equations, RAIRO Anal. Numer. 13 (1979), no. 3, 

201-226. MR MR543933 (81c:65044) 
[BJR90] Alain Bamberger, Patrick Joly, and Jean E. Roberts, Second-order absorbing boundary 

conditions for the wave equation: a solution for the corner problem, SIAM J. Numer. 

Anal. 27 (1990), no. 2, 323-352. MR MR1043609 (91b:35066) 
[BJT00] E. Bcache, P. Joly, and C. Tsogka, An analysis of new mixed finite elements for the 

approximation of wave propagation problems, SIAM J. Numer. Anal. 37 (2000), no. 4, 

1053-1084 (electronic). MR MR1756415 (2001d:65124) 
[BR99] W. Bangerth and R. Rannacher, Finite element approximation of the acoustic wave 

equation: error control and mesh adaptation, East- West J. Numer. Math. 7 (1999), 

no. 4, 263-282. MR MR1738435 (2000i:65148) 
[BR01] Wolfgang Bangerth and Rolf Rannacher, Adaptive finite element techniques for the 

acoustic wave equation, J. Comput. Acoust. 9 (2001), no. 2, 575-591. MR MR1853643 

(2002f: 76049) 

[BS01] Ivo Babuska and Theofanis Strouboulis, The finite element method and its reliabil- 
ity, Numerical Mathematics and Scientific Computation, The Clarendon Press Oxford 
University Press, New York, 2001. MR MR1857191 (2002k:65001) 



L°°(L 2 )-NORM A POSTERIORI BOUNDS FOR WAVE EQUATION 



17 



[BS05] Christine Bcrnardi and Endrc Sli, Time and space adaptivity for the second-order wave 
equation, Math. Models Methods Appl. Sci. 15 (2005), no. 2, 199-225. MR MR2119997 
(2005j:65105) 

[BV04] Christine Bcrnardi and Rdiger Verfrth, A posteriori error analysis of the fully dis- 
cretized time- dependent Stokes equations, M2AN Math. Model. Numer. Anal. 38 
(2004), no. 3, 437-455. MR MR2075754 (2005g:65131) 

[CB02] Carsten Carstensen and Sren Bartels, Each averaging technique yields reliable a pos- 
teriori error control in FEM on unstructured grids. I. Low order conforming, non- 
conforming, and mixed FEM, Math. Comp. 71 (2002), no. 239, 945-969 (electronic). 
MR MR1898741 (2003e:65212) 

[Cia78] Philippe G. Ciarlet, The finite element method for elliptic problems, North-Holland 
Publishing Co., Amsterdam, 1978, Studies in Mathematics and its Applications, Vol. 
4. MR 58 #25001 

[CJT93] Gary Cohen, Patrick Joly, and Nathalie Tordjman, Construction and analysis of higher 
order finite elements with mass lumping for the wave equation, Second International 
Conference on Mathematical and Numerical Aspects of Wave Propagation (Newark, 
DE, 1993), SIAM, Philadelphia, PA, 1993, pp. 152-160. MR MR1227833 (94d:65058) 

[CKNS08] J. Manuel Cascon, Christian Krcuzer, Ricardo H. Nochetto, and Kunibcrt G. Sicbcrt, 
Quasi-optimal convergence rate for an adaptive finite element method, SIAM J. Numer. 
Anal. 46 (2008), no. 5, 2524-2550. MR MR2421046 (2009h:65174) 

[DS81] Vassilios A. Dougalis and Steven M. Serbin, On the efficiency of some fully dis- 
crete Galerkin methods for second-order hyperbolic equations, Comput. Math. Appl. 7 
(1981), no. 3, 261-279. MR MR614183 (82e:65106) 

[D96] Willy Drfler, A convergent adaptive algorithm for Poisson's equation, SIAM J. Numer. 
Anal. 33 (1996), no. 3, 1106-1124. MR MR1393904 (97c:65139) 

[EJ95] Kenneth Eriksson and Claes Johnson, Adaptive finite element methods for parabolic 
problems. II. Optimal error estimates in Loo£2 o-nd Loo^oo, SIAM J. Numer. Anal. 
32 (1995), no. 3, 706-740. MR 96c:65162 

[GLV08] Emmanuil H. Georgoulis, Omar Lakkis, and Juha M. Virtancn, A posteriori error 
control for discontinuous Galerkin methods for parabolic problems, preprint submitted 
to Journal, under revision 0804.4262v2, arXiv.org, 2008. 

[HS01] Paul Houston and Endrc Sli, Adaptive Lagrange- Galerkin methods for un- 
steady convection- diffusion problems, Math. Comp. 70 (2001), no. 233, 77-106. 
MR MR1681108 (2001f:65114) 

[Joh93] Claes Johnson, Discontinuous Galerkin finite element methods for second order hyper- 
bolic problems, Comput. Methods Appl. Mcch. Engrg. 107 (1993), no. 1-2, 117-129. 
MR MR1241479 (95c:65154) 

[KM05] Ohanncs Karakashian and Charalambos Makridakis, Convergence of a continuous 
Galerkin method with mesh modification for nonlinear wave equations, Math. Comp. 
74 (2005), no. 249, 85-102 (electronic). MR MR2085403 (2005g:65147) 

[LM06] Omar Lakkis and Charalambos Makridakis, Elliptic reconstruction and a posteriori 
error estimates for fully discrete linear parabolic problems, Math. Comp. 75 (2006), 
no. 256, 1627-1658 (electronic). MR MR2240628 

[Mak92] Ch. G. Makridakis, On mixed finite element methods for linear elastodynamics, Numer. 
Math. 61 (1992), no. 2, 235-260. MR MR1147578 (92j:65142) 

[MN03] Charalambos Makridakis and Ricardo H. Nochetto, Elliptic reconstruction and a pos- 
teriori error estimates for parabolic problems, SIAM J. Numer. Anal. 41 (2003), no. 4, 
1585-1594 (electronic). MR MR2034895 (2004k:65157) 

[MN06] , A posteriori error analysis for higher order dissipative methods for evolution 

problems, Numer. Math. 104 (2006), no. 4, 489-514. MR MR2249675 (2008b:65114) 

[Pic98] Marco Picasso, Adaptive finite elements for a linear parabolic problem, Comput. Meth- 
ods Appl. Mech. Engrg. 167 (1998), no. 3-4, 223-237. MR MR1673951 (2000b:65188) 

[Stc07] Rob Stevenson, Optimality of a standard adaptive finite element method, Found. Com- 
put. Math. 7 (2007), no. 2, 245-269. MR MR2324418 (2008i:65272) 

[S96] Endrc Sli, A posteriori error analysis and global error control for adaptive finite 

volume approximations of hyperbolic problems, Numerical analysis 1995 (Dundee, 
1995), Pitman Res. Notes Math. Scr., vol. 344, Longman, Harlow, 1996, pp. 169-190. 
MR MR1405623 (97d:65057) 



18 



E.H. GEORGOULIS, O. LAKKIS, AND C. MAKRIDAKIS 



[S99] , A posteriori error analysis and adaptivity for finite element approximations of 

hyperbolic problems, An introduction to recent developments in theory and numerics 
for conservation laws (Freiburg/Littenweiler. 1997), Lect. Notes Comput. Sci. Eng.. 
vol. 5, Springer, Berlin, 1999, pp. 123-194. MR MR1731617 (2001d:65119) 

[Ver96] Rdiger Verfrth, A review of a, posteriori error estimation and adaptive mesh-refinement 
techniques, Wiley- Teubner, Chichester-Stuttgart, 1996. 

[Ver03] R. Verfrth, A posteriori error estimates for finite element discretizations of the heat 
equation, Calcolo 40 (2003), no. 3, 195-212. MR MR2025602 (2005f:65131) 

Department of Mathematics, University of Leicester, University Road, Leicester 
LEI 7RH, United Kingdom. 

E-mail address: Enraianuil.GeorgoulisOnics.le.ac.uk 

Department of Mathematics, University of Sussex, Falmer near Brighton, GB-BN1 9RF, 
England. |http: / /www. maths. Sussex. ac.uk/Staff/OL| 
E-mail address: 0.Lakkis@sussex.ac.uk 

Department of Applied Mathematics, University of Crete, GR-71409 Heraklion, Greece; 
and Institute for Applied and Computational Mathematics, Foundation for Research 
and Technology-Hellas, Vasilika Vouton P.O. Box 1527, GR-71110 Heraklion, Greece. 

E-mail address: makratem.uoc.gr 



